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ABSTRACT 

subsystTiuh Carbon ' sond°am?ne water desorbed 
subsystem called SAWD. 

This manual describes ^■« 1 “; 8 :l!i ,1 ^ti"d.«Ti!S: and/or" each 
Emul ation/Simul ati on Computer .Hodel- P°J * nd chem(ca l laws which govern 

individual component, the fundament _ physi^ ^ stated , an d when necessary, 

data'i s presented support empirically developed relationships. 
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1.0 

INTRODUCTION 

1.1 

Background 


ESCM is a computer model which was developed to demonstrate the 
utility of a major portion of the Emulation/Simulation Sizing, and 
Technology Assessment Program (ESSfAP). See Reference 1. ESSTAR- 1 s 
a concept for software tools that will support the total engineering 
process of the Space Station beginning with concept definition and 
continuing on through mission operations. The interaction of the 
software tools with the design and operational phases of a flight 
system is shown in Figure 1. The philosophy of ^ software 
concept is to have the analysis software for each step in the design 
process precede development of hardware in order to provide the 
greatest design and cost benefits. 

Of the many systems in the Space Station, the Environmental Contro 1 
and Life Support System (ECLSS) was selected to demonstrate the 
utility of the ESSTAP concept because of its complexity (number of 
components and dynamic operational capability), aval 1 abi 1 ity o 
operational data for checkout, and growth potential. The ESC 
program was targeted to evaluate and demonstrate the benefits o 
phases 2 through 4 software. For this purpose, an evaluation of the 
six major ECLSS functions was performed as published in the program 
document "ESCM-EM-02" (Reference 2). As a result of this study, a 
subgroup of the Air Revitalization Subsystem (ARS) shown in Figure 2 
was selected for modelling because of its dynamic complexity, growth 
potential, design tool utility, and the independent nature of the 
ARS compared to other subsystems. 

The resulting ESCM model has been used to demonstrate the use of a 
computer program in the design, development, and test of applicable 
development hardware. Specifically, through the use of the program, 
design hardware and software can be verified. Failure Mode and 
Effect Analysis (FMEA ) can be assisted, and test planning can be 
technical ly improved. 


1.2 


Overview of Modelled System 


fhe ARS functions and components are depicted in Figure 2. The 
system to be simulated involves the Space Station cabin, crew, a 
heat exchanger for temperature control , a heat exchanger for 
humidity control, fans, a SAWD (Solid Amine Water Desorbed) carbon 
dioxide removal system, and a two-gas controller for oxygen 
replenishment and total pressure control. 


The cabin is a volume which houses the crew, two-gas 
certain vehicle heat loads representing equipment, 
of the cabin is also modelled. The crew breathes in 
cabin air and releases sensible heat, latent heat 
respiratory water vapor, and carbon dioxide. 


controller, and 
Air leakage out 
oxygen from the 
in the form of 
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Figure 1 


Application Of Software 
Operational Phases 


Tools To The Design 
Of A Flight System 


And 
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Figure 2 

ARS - SAWD II Schematic For Demonstration Model 
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nvprview of Modelled System (Continued) 


Cabin air is drawn by fans g i ntc . three 

passes the air through a heat exc™ 1 ^ ff its heat to 

from the. air. In the heat exchanger . tt. of 60«F. The 

water which enters the he controls the flow of water to 

" • 

isrjs jus •; s-H rr«r7«£S t 

to water which enters the h * a * . e £? h Ju? s exchanger, and a portion of 

vapor in the air is condensed in this excnang^, ^ ^ a wate r 

the main air stream, with c ° ndens ® . P removed. This humidity 

separator wher % th ® ^o^cei ves air flow directly from the SAW D 
control heat exchange! r , al e s ° t p® Thus the effect of any residual 

moisture i0 from ^SAWD* Wem oh ’cabin relative humidity is 

lessened. 

.. hin for use by the SAWD carbon 
The third fan draws ai r from removal system contains two 

dioxide (CO 2 ) removal syste . accumulator, a water pump, a 

amine beds, solenoid valves % and instru- 

steam generator, a c °" tr ° system operates to remove CO 2 contin- 
mentation. In general, ^ . des0 rbing each bed. During 

uously by alternately abs 9 . oar ticles of the bed. During 

’ thr^h 

amine” 6 . n d* «rVy tV ^ to the accumulator. Figure 3 shows a 

cross section of the bed canister. 

A,r from the SAWD bed “"^olos “SeT The^air 

control heat exchanger where any moist ^ Cabln 

returning to the cabin from the SAWD system row ^ contro))ed by a 

air total pressureandorygenpat J J 4>ns 0(ygen as required to 

Shuttle type two ^s contro, e 3 . 0 9 tQ 3 <23 psia . It 

maintain the oxygen partial pressu tQ 14<813 ps ia only 

admits nitrogen to bring ^‘J^’J^een brought above 3.23 psia. 
after the oxygen partial pressure 


MODELLING OF ARS SUBGROUP FOR G189A 


- ■ 

The ARS subgroup to be analyzed is ^ ow J o in b g 1 ^ U o r d e 1 2 1 * ed i n terms of 

s ub rout i ne s^'and ' requirement s^ oT ^the G189A computer program which is 
the basis for the ESCM computer program. 


4 
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modfi LING OF APS SUBGROUP FOR G189A (Continued) 

The modelled system shown in Figure 2 , and heat 

following component models are 

needed : 


(7) Pump 

(8) Steam Generator 

(9) Splitter 

(10) Mixer 

(11) All Purpose Component 


(1) Crew 

(2) Cabin 

(3) Fan 

(4) Heat Exchanger 

(5) I R 45 Bed 

(6) Tank 

By combining these components i n Figure 2. 

Figure 4 is generated to model the AKb SU y the H ESCM item number and 
Table 1 presents a prison of these two figures 

the corresponding subroutine. Qf all> the crew which 

reveals some subtle dlffe ^ en carbon dioxide and water vapor is 
consumes oxygen and produ 1 } that contains the sensible 

placed in series at the en tQ keep the schematic simple, 

heat exchanger and fan* . Eo D i ace the crew in parallel with 
An alternative would have been but this setup would have 

the sensible heat exchanger ^ nd ’ d the reS ult to the cabin 

required an extra splitter an ’ f car b 0 n dioxide and water 

fan fl ow. 

Another difference V c'jr'lhe G18W J«puUr° 

the secondary. For the cabin , the P J crew . The sec0 ndary 
contains the sensible heat ex g and the humidity 

flow path, therefore, must contain^both thej ^ u$ing a sp1itter 

control heat exchanger. f the P fi 0 w to the SAWD and the 

(Item 3) to direct a portion of t . . The calculation 

remainder to the hum l^ 1 1 t ^ c d 5 S c USS ed in Section 3.0. The return of 
of the split ratio wl1 J be .. t the humidity control heat 

the SAWD flow enters the inlet line to 

exchanger in both figures. 

Another difference occurs ' " * he , "“the ' cart'on^d i oxi de removal 
crosses, and tees which «’ st J" t V E SCH by splitters or 

system. Tees in 1 ines felled in the^tbtn y ^P and then 

mixers depending on whether a pnterinq a tee and then mixing, 
splitting or whether two flowsare( ( g sp1itte r) entering the 

Accordingly, item 14 represents ; the , tee U sp ^ ^ (a mixer) 

™ :^ve selects one 






Figure 4 

ARS - SAWD II Schematic As Implemented In 
ESCM Computer Program 
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Table 1 

fscm Component Number - Subroutine Cross Reference Tab le 


Component Description 

Crew in cabin 
Cabin 

Split to humidity control HX or SAWD 
Cabin sensible HX fan 
Cabin sensible HX 

Cabin sensible HX cooling fluid boundary 
Cabin condensing HX 
Condensing HX cooling fluid boundary 
Split to water separator or bypass 
Water separator 

Mix water separator and bypass 

Condensing HX fan 

Mix SAWD and condensing HX flows 

Split to SAWD bed #1 or #2 

Mix steam or cabin air to bed #1 

SAWD bed #1 

Split to cond. HX or C0 2 recovery 
Mix bed #1 and bed #2 flows 
SAWD fan 

Mix steam or cabin air to bed #2 
SAWD bed #2 

Split to cond. HX or C0 2 recovery 
Mix bed #1 and bed #2 exit gases 
Split to preheat or C0 2 accumulator 
Mix steam generator and preheat streams 
Split to SAWD bed #1 or #2 
Water supply tank 
Water supply pump 
Steam generator 

Split desorbed gas to cabin or C0 2 accui 
Mix cabin air and SAWD air 
Carbon dioxide accumulator tank 


Subroutine 

Name 

Subroutine 

No. 

SUITS 

2 

CABIN 

1 

SPLIT 

10 

FAN 

23 

ANYHX 

4 

ALTCOM 

49 

ANYHX 

4 

ALTCOM 

49 

SPLIT 

10 

ALTCOM 

49 

GASMIX 

6 

FAN 

23 

GASMIX 

6 

SPLIT 

10 

GASMIX 

6 

IR45 

73 

SPLIT 

10 

GASMIX 

6 

FAN 

23 

GASMIX 

6 

I R 45 

73 

SPLIT 

10 

GASMIX 

6 

SPLIT 

10 

GASMIX 

6 

SPLIT 

10 

TANKG 

30 

PUMP 

22 

SMGEN 

27 

,r SPLIT 

10 

GASMIX 

6 

TANKG 

30 
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Mnnpl l inr. OF ARS SUBGROUP FOR G189A. (Continued) 

of two entering flows to be passed ^ough^the J^are "model 1 ed as 

the three-way valves upstrea , three-way valves downstream of 
mixers (items 15 and 20); while the three way^va ^ and 22) The 

each amine bed are model 1ed h iLen the two amine beds in the Figure 
cross, valve, and two tee s bet mixers an d splitters 
2 schematic are modelled as a co b i t mix tw0 flows together and 
The action of a cr ? ss .° *» to one of two distinct paths, can be 
ri^rr tToVlf a mixer and a setter. 

The last difference occurs in e ** 

recovery. The bac n .^t^v^ wh?^ regulate the pressure in the 
regulators, onfl £ e ’ a t 0 f t he accumulator are not modelled as 

accumulator and the flow oub ° f , th ® h ^ action is input in the form 

separate components; but 1 "1* ' fl ’ 0 * leaving the item 32 accumulator, 
of equations which compute the flow leaviny 

modelling of components 

As discussed in Section 2.0, o y subqroup. The following 

need be available to a . nal ^ ze a i h ytica i models, the generation of 
sections will discuss these analyt^ papaiI 4 t:ers necessary to 

SKe ^"al. model In 

^ofe cases^the^ reader 1 will ^e preferred to the G189A manual for a 
description of the analytical model. 

i. thp control of cabin temper- 
In addition to the ® leven cabl ! n oxygen partial pressure, and 

ature, cabin air total P^^ r .’ gpoLYl and GP0LY2 subroutines. A 

description ^fT he^ SffrM be provided. 

Crew 

The crew affects the gas strain tempe ^ at ^VesulTof^its metabolic 

dioxide content, and water content as a result^ ^ tQtal 

functions. To solve for these influen^ shown ip Figure 5 

metabolic load schedule i d different loadings for eating, 
shows a probable profile . with dl ”f e average 1oading over the 24 
working, exercising, and s P 9* equivalent to a carbon 

hour period Is 477 Btu/man-hour This ‘ Fronl the crewman 

di0 tabollc 9 To e ad a !'?!;e per man sensible' and latent splits are calcu- 
lated In GP0LY1 using the following relations. 


0, - 430 * (HR/1000 + 10) * ( TCAB 

Suilll - 0.22 * MR t *.* * («« - «•> 

QL = Larger of QL and QLMIN 

QS = MR - QL 


60) 



METABOLIC RATE PER MAN (BTU/HR) 
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Figure 5 

Crewman Metabolic Heat Load Profile 
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3.1 Crew (Continued) 

Where: MR 

TCAB 

QL 

QLMIN 

QS 
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3.2 


per man metabolic heat Joad, Btu/man-hr 
cabin air temperature, °F 
latent heat added to the gas stream 
per man, Btu/man-hr 
minimum latent heat added to the 
qas stream per man, Btu/man-hr 
sensible heat added to the gas stream 
per man, Btu/man-hr 


With the sensible and latent heat loads known, subroutine SUITS is 
used to compute the following: 


0T 

WC02 

W02 

WH20 


total gas stream heat load from the crew, Btu/hr 
total crew carbon dioxide generation rate, lbm/hr 
crew oxygen gas usage, lbm/hr 
crew water vapor generation rate, lbm/hr 


which is not described in Reference 3: 


T e = 

Where 


(m n - c p Ti + N * QS) / (mi c p ) 


'e 

mi 


Ti 

N 

QS 


exiting temperature of f low T° wh l£ h 
the crew sensible heat is added, F 
flow to which crew sensible heat is 
added, lbm/hr 

specific heat of flow, Btu/lbm- F 
inlet temperature of flow, °F 
number of people in crew 
sensible heat added to the gas 
stream per man, Btu/man-hr 


iss’« 

against design criteria and (2) if steady state appropr1ate 

analyses are to be performed. TJe »pt ions selected as ppp 
for ESCM are to perform steady state analyses ana not to H 
comparisons against design criteria. 

Cabin 

I • il 0 prrM -i c essentially a( mixing chamber for 
It ?? b ” H the^sensible^eat J, changer, condensing 

heat exchanger, the crew, the SAWD cart d l»" de s ^° r v a ^ 
and the oxygen and nitrogen supply systems. J 

ttste Ssed for steady state and transien t calcul 

steady state performance of the ^ utioos are performed in 

K analyses performed 

in these subroutines is given in Reference 3. 


11 
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Cabin (Continued) 

The following options are used in the ESCM application: 

(1) The primary scarce flow pressure is used as the cabin mixture 
pressure during steady state calculations. 

(2) Reset the cabin primary and secondary side outflow data equal 
to the initial input data for five steady state passes. 

(3) Calculated entrained liquid in the cabin will remain in the 
cabin air mixture. 

(4) Cabin outflows will be calculated for both the primary and 
secondary paths. 

(5) Mass addition temperatures are specified. 

and nitrogen addition rates input by the two gas controller. 

(7) The inner wall is considered to be a boundary temperature. 

(8) The furnishings are considered to be a boundary temperature. 

(9) The outer wall is considered to be a boundary temperature. 
Input specific to the ESCM application are as follows: 

(1) Initial flows, temperatures, and pressures: 


Cabin gas exit temperature, F 
Inlet pressure, psia 
Exit pressure, psia 
Water vapor flow, lbm/hr 
Oxygen flow, lbm/hr 
Nitrogen flow, lbm/hr 
Carbon dioxide flow, lbm/hr 


Primary 

70 

14.7 

14.7 

89.7 
2046.8 

7199.4 

63.4 


Secondary 

70 

14.7 

14.7 

13.0 

303.76 

1068.67 

9.45 


Control setpoints, tolerances, and limits: 

Cabin gas temperature control setpoint, °F 
Cabin gas temperature control tolerance, F 
Cabin gas total pressure control setpoint, psia 
Cabin gas total pressure control tolerance, psia 
Cabin gas oxygen pressure control setpoint, psia 
Cabin gas oxygen pressure control tolerance, psia 
Cabin gas dew point control setpoint, F 
Cabin gas dew point control tolerance, F 
Maximum allowable CO 2 pressure, mmHg 

Maximum allowable trace contaminant level , ppm 

•These ve.ee, .he -ere., “ 

setpoints and tolerances wh ch are hardwired 3>12>u 

Accordingly, "a ch.Si " these inputs has no influence on the two 
gas control setpoints and tolerances. 


12 
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Cabin (Continued) 

(3) Mass addition rates and temperatures. 


Water vapor addition rate, lbm/hr 
Entrained water addition rate, lbm/hr 
Water vapor addition temperature, F 
Oxygen addition temperature, F 
Nitrogen addition temperature, F 
Carbon dioxide addition temperature. 

Trace contaminants addition temperature. 
Special flow no. 1 addition temperature, F 

0.233 

0.000 

70 

70 

70 

70 

70 

70 

Net flows into cabin resulting from primary and 
flow stream conditions: 

secondary 

Net carbon dioxide addition, lbm/hr 
Net trace contaminants addition, Ibm/nr 
Net special flow no. 1 addition, lbm/hr 

0.0 

0.0 

0.0 

Miscellaneous cabin parameters. 


Equipment, lighting, and miscellaneous 
heat loads, Btu/hr 
Outboard leakage rate, lbm/hr 
Cabin free volume, ft 

17065 

0.0833 

8000 


The above listed values are determined as follows 
(1) Water vapor addition rate: 


From Reference 2 Table 1, 
Hygiene latent load 
Food preparation latent load 
Experiments latent load 
Laundry latent load 


=0.96 lbm/man-day 
=0.06 lbm/man-day 
=1.00 lbm/day -module 
= 0.13 lbm/man-day 


Assume number of men - 


Hygiene 

Food 

Experiments 

Laundry 


= 4 x 0.96 
= 4 x 0.06 
= 1.00 
= 4 x 0.13 


4 


for conservative sizing 


3.84 

0.24 

1.00 

0.52 

5.60 lbm/day = 


0.233 lbm/hr 


(2) Equipment, lighting, and miscellaneous heat loads: 

Fran Reference 2 Table 1, the Non-ECLSS thermal load is 5.00 kw 
per module. Therefore: 

Q = 5.00 x 3.413 = 17065 Btu/hr-modul e 
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Cabin (Continued) 

(3) Outboard leakage rate: 

_ o T , hlo i the cabin atmosphere leakage is 2.00 

^day-Se- Accordingly , the per hour leakage rate is 
0.0833 1 bm/hr-modul e. 


( 4 ) 


( 5 ) 


'abin free volume: 

*ny volume can be input by the user for the cabin volume. An 
assumed volume of 8000 cubic feet is used. 

Control setpoints are taken as roughly p osi t ions°specified in 
for allowable temperatures and gas composition* v 

Roforpnce 2. 


Ref. 2 
Range 


Sel ected 
Value 


Cabin gas total pressure, psia 
Oxygen pressure, psia 
Water vapor pressure, psia 
Dew point, °F 

Carbon dioxide pressure, mmHg 
Dry bulb temperature, °F 


10.0-14.7 
2. 7-3. 2 
0.12-0.26 
40-60 

3 avg 7.6 max 
65-75 


14.7 

2.9 

0.20 

50 

7.6 max 
70 


(6) Initial flows: 

^5 SSte“?.°ibSTnS b S tWlSSn, 

calculations: 


.02 


= 14.7 psia 
= 2.9 psia 


N2 


= 14.7 - 2.9 - 0.2 - 0.06 
= 11.54 psia 


H20 = 0.2 psia 


P C02 = 3 '° mmH9 


0.06 psia 

F ;: section 3.3, the sensible and condensing heat exchanger 
fan flows are: 


Sensible heat exchanger 
Condensing heat exchanger 


2100 

300 

2400 cfm 


14 
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3.2 Cabi n (Continued) 


Gas 

Mw 

Press 

(psia) 

Density* 

(PCQ 

Flow, 

Primary 

lbm/hr 

Secondary 

Ni trogen 

28.016 

11.54 

.0569 

7169 

1024 

Oxygen 

32 

2.90 

.0163 

2054 

293 

Water vapor 

18.016 

0.20 

.000634 

80 

11 

Carbon dioxide 

44.01 

0.06 

.000449 

57 

8 

144PM, 
* P = TT 

- where T = 70 

°F = 530°R 




The above values were used as the first guess. After having 
run the ESCM program, the program computes the actual steady 
state flows. These computed values can then be used as a 
better first guess for subsequent cases. These better first 
guesses are the values shown in Reference 4 Table 8 for the 
initial flows in the input specific to ESCM. 


3.3 Fan 

Three fans exist in the ARS subgroup modelled by ESCM. They are the 
sensible heat exchanger fan, condensing heat exchanger fan, and the 
SAWD carbon dioxide removal system fan. The G189A subroutine FAN is 
used to model the fans. Reference 3 gives a complete description of 
the analytical models used in the subroutine. Of the different 
analytical models available in that subroutine, the following is 
used for the sensible and condensing heat exchanger fans: 

(m T ) = 60 * cfm * p 

Pe = p i 

Q = w * 3.413 

T - t. + 0 

e 1 (>y e c p 

R = (m T ) e /(m T ) i 

Kc>e = (m nc^i * R 

(\) e = (" v )i * R 

(m ] ) e = (n» 1 ) i * R 

( m 02 ) e = ( m 02>1 * R 
H^e = (m N2 } e * R 
^ m C02^e = ^ m C02^e * R 


15 


SVHSER 9504 



TECHNOLOGIES 


3.3 


Fan (Continued) 
Where: 


m = mass flow, lbm/hr 

r : ^ 

l : 1 If 


3.3.1 


Subscripts: 


j = total 
e = outlet 
i = inlet 

nc = non-condensables 

v = water vapor 

1 = entrained liquid 

02 = oxygen 

N2 = nitrogen 

C02 * carbon dioxide 


The only input needed is 


the fan cfm and power 


for each fan. 


rrmH^nsing Heat Exchanger Faji 


For the condensing 
used to obtain the 


heat exchanger fan, the 
fan cfm and power: 


following 


analysis was 


Assume : 

dry bulb temperature, T^b 70 F 
relative humidity, RH = 55% 


Therefore: 

o A 70°F = 0.3631 psia 

P =1 70°F 55% RH = 0.55 . 0.3631 - 0.1997 ps. a 

The corresponding dew point is obtained from the stean 
Reference 5. 


Tdp - Tsat * °- 1997 P $1a * M - 1 ° F 


OJ 


Mv Pv 
M a p a 


My p v 

M a p T- p v 
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teSSologies 


Condensing Heat Exchanger Fan (Continued) 

Where : 

« = specific or absolute humidity, ]bm vapor/lbm air 

M v = molecular weight of the Jbm/lb-mole 

M l = molecular weight of air, lbm/lb-mole 
P 3 = absolute pressure of vapor, psia 
Pg = absolute pressure of air, psia 

Therefore: 

18.016 * _0a997__ , o. 01 268 
i = 28795" 10-0.1997 

In the condensing heat exchanger, assume: 

(T H 20)i = 45 ° F 

(T a )i - 70°F 

Effectiveness, = 0.85 

(T a ) = 70 - 0.85 (70 - 45) = 48.8°F 

Assume saturated conditions at exit. Therefore from Ref. 5: 

T sat = 48.8°F P sa t = ( p H20) e = °- 1703 P sia 
18 * 0.1703 „ 0.01078 


w e = 28.92 * (10-0.1/03) 


= W1 .w e = 0.01268 - 0.01078 = 0.00190 

The mass flow of dry air is determined from. 

_ m l 
m a wi n - w ou t 

Where: 

m a = mass flow of dry air, lbm/hr 

mi = mass flow of condensed water, lbm/hr 

The mass flow of water to be condensed is determined for a four man 
crew. From Reference 4 Table 1: 


Load 

Sweat and Respi ration 
Hygiene Latent Load 
Food Preparation 
Experiments 
Laundry 

Average SAWD Load 
Total 


= 8 * 4.02 = 32.16 Ibm/day 

= 8 * 0.96 = 7.68 lbm/day 

= 8 * 0.06 = 0.48 lbm/day 

=1.00 = 1-00 lbm/day 

=8 * 0.13 = 1.04 lbm/day 

= 8 * 4.88 = 39.04 lbm/day 

81.40 lbm/day 
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3.3.1 


Condensing Heat Exchanger Fan (Continued) 


For a 4 man system 
water condensation 
of dry air yields: 


the load is 40.7 lbm/day. 
rate into the above equation 


Substituting this 
for the mass flow 


"a -gmm * 892 - 5 ' bl " /hr 


The fan cfm then is: 

, m a 
cfm “ JTO 

Where: 


cfm = fan flow, cfm 


f 

cfm 


density of air, lbm/ft 3 

5 , = 292 cfm Use 300 cfm 
.051 * bu 


To compute the fan power, the following equation and assump 
tions are used: 


cfm *AP 

w = 8.5 


Where: 

w = fan power, watts 

AP = fan head rise, In-HgO 

tj = fan efficiency 

cfm = fan volumetric flow, cfm 


Assume : 

AP = 2.0 In-H20 
T = 0.5 

292 * 2.0 _ •) ^7 watts 

w = 8.5 * O ■ 137 wa 

The condensing heat exchanger sensible heat load is: 
q = m c p AT = 892.5 x 0.24 x (70 - 48.8) = 4541 Btu/hr = 1330 watts 
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3.3.2 Sensible Heat Exchang er Fan 

The loads for the sensible heat exchanger are as follows: 

Non-ECLSS loads = 5000 w/module 
Modelled ECLSS loads: 


CHX Fan 

SAWD Fan + Steam 


= 137 watts 
= 86 watts 

223 watts 


7000 * 4 = w 

For 4 man crew = -3.413 * 2"4 

£w - 6000 + 223 + 342 - 1330 - 4235 watts 

The value of 1330 watts is subtracted^in ^Tontoslng" 
heat C exchanger. ^Fo^the sensible heat exchanger assume the 
foil owing : 

(T a ). = 71 ,8°F (1.8°F additional for fan) 

(T H 20)i = 60 ° F 

effectiveness, < = 0.80 

(T a ) e - da), - Ohzo),] 

= 71.8 - 0.80 (71.8 - 60) 

= 62.4°F 

The flow of air is then given by. 

Q 423 5 * 3.413 


c^ST = 0724 - (71 .8 - 6274) 


= 6407 Ibm/hr 


cfm 


"a _ 


6499 


,60 ’ .051 * 60 


= 2093 Use 2100 cfm 
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3.3.2 


3.3.3 


Sensible Heat Exchanger Fan (Continued) 
The fan power required is: 


w 


cfm * AP 
8.5 * t? 


Assume: 

AP = 2 In-H 2 0 

-q =0.60 

2093 * 2.0 
w = 8.5 * O.bU 


= 821 watts 


SAWD Fan 

For the SAWD fan, the following analytical model is used: 
(m T ) e = ( m j)'j 


p e 

= Pi 

• 

Q 

= w * 3.413 

T. + 9 

^e 

= Tl + ^T~ 

(m ) 
v ne'e 

= Kc } i 

("\,) e 

= K>i 

( m l)e 

= ( m l)i 

( m 02^e 

= ( m 02 ) i 

( m N2^e 

= ( m N2^i 

( m C02^e 

= ( m C02^ i 

The SAWD fan 

flow however 


and then into the second bed which is dry. 
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3.3.3 


3.4 


SAWD Fan (Continued) 

The fan cfm is shown in Figure 6 as a function of time within a 
phase. During energy transfer where air is flowing through two 
beds, the fan flow follows the energy transfer phase curve shown 
in Figure 6. When the energy transfer phase is completed, the 
flow which existed at the end of the energy transfer phase is 
used as the starting point flow on the absorption phase curve 
and a pseudo-el apsed time is created. The fan cfm then increases 
along the absorption phase curve per this pseudo elapsed time. 

For the SAWD fan power, assume: 

cfm =21 
AP = 7 In-H 2 0 
i? = 0.35 

p 0wer = Q ..AP = flJl-Z— = 49.4 watts Use 50 watts 
rower 8>5t] 8.5 * 0.35 


Heat Exchangers 

Two heat exchangers are in the ARS subgroup modelled by ESCM. They 
are the cabin air sensible heat exchanger and the cabin air humidity 
control heat exchanger. The subroutine used in ESCM to model the 
heat exchanger is called ANYHX and a complete description of all the 
equations and options used in the subroutine is given in Reference 3. 
The option used in the ESCM is to input the heat exchanger UA and 
then to use the standard steady state effectiveness - NTU relations 
for a counter flow heat exchanger to compute the heat exchanger per- 
formance and exit temperatures. 


The UA for the sensible heat exchanger is 6240 Btu/hr-°F and is 
calculated as follows: 


(^H20) O ut m Cp + ^H20)-j n 


(4235 + 821) * 3.413 + 6Q 
1900 


= 69. 1°F 


c min 

c max 


6500 * 0.24 
‘ 1900 * 1.0 


0.821 


For a counter flow heat exchanger with € 
0.821, NTU = 4.0. 


0.80 and c m -j n /c max - 


Therefore : 

UA = NTU * c min = 4.0 * 6500 * 0.24 = 6240 Btu/hr-°F 
This is a constant UA for the the sensible heat exchanger. 
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Figure 6 
SAWD Fan Flow 
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3.4 


Hpat Exchangers (Continued) 


The t UA ne fO GP0m 

routine GP (0.52 + 0.00026 * mn20) 


UA = e 


(2.809 - .00169 * m H 20) * m ma x 


3.5 


3.5.1 


Where : 

. ay is the maximum of the following: 

(1) ma1r . total air flow into heat exchanger including the 
1 w vapor flow. 


(cp)eff 

(2) m ef f = "air '(c p )ai7 


The effective specific heat is defined in Reference 3 in the 
subroutine description for ANYHX. 

SAWD Bed . 

The analysis of the.SAWO b«M ^"stb^ineU^ly 

the name of the resin used in the : bed. Refere „ce 3 6189ft 

written for ESCH and is "°‘ ^described the analysis modelled in 

subroutine* 1 IRA^t^compute 1 the carbon dioxide absorption and desorp- 

tion from the resin. 

As shown in Figure 2, 

header volumes and a centra . this central volume is 

with dry amine. In the computer -model .Uns s « thfi model 

divided into five segments. Th , t and exit header, five void 

of the SAWD canister consists or an f oa ni pluS can 

volume segments and five ^^nsor The air flows into the inlet 
segments, and a temperature one through five, out the outlet 

header, through void er ature 9 sensor. In the void seg- 

header, and lastly past the tempera d with the corres- 

ments, water and carbon sec ?ions describe the equa- 

ponding resin segments. occurring in the headers, the b 

tions which model the process^ sensor, 

void and resin segments, and the tempera 

The Headers 

The fundamental process m0 ^ 1 p 1 e ^ 1 "esi di ng* i ^the 5 header!^ This 
incoming gas with the gas J libr ium temperature and a requ' ire 

mixing process results in an eq Dressur e less than or equal to 

amount of exiting gas to ma^nta a P mixture of non-condensables 

the downstream. If the mxeu y« 

and vapor, the following analysis is used. 
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SAWD canister is 


divided into seven segments 


including 5 bed segments 


I 

n_ 

I 

e 

t 


e 

a 


[Segment 

#1 

void 

volume 


|Segment 

#2 

void 

volume 


Segment 

#3 

void 

volume 


Segment 

#4 

void 

volume 


Segment 

#5 

void 

volume 



[Segment 

#1 

resin + 
foam + 
can 


Segment! 



[Segment 

1 #3 

resin + 
foam + 
can 



0 

— i u 

JTemp^t 

Isen- | 
sor p 


Figure 7 

SAWD Canister Simulation 
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The Headers (Continued) 

HA ■ ( M „c c p,"C , i T ’ + (M v>i "S’* + (M,) < 

HS ■ ( M nc c p.nc>ST T ST + (M *>ST h g,ST + (M 1 } ST h l ,ST 


HMIX = HA + HS 


Where : 


Total enthalpy of gas entering header, Btu 

Total enthalpy of gas already stored in header, Btu 

Total enthalpies of gases entering and stored, Btu 


Subroutine HBALNC is called to establish the mature temper- 
ature mass of vapor, and mass of entrained liquid for an 
initial starting pressure equal to the downstream Pr^sure. 
With this new mixture temperature, mass of vapor, and mass or 
entrained liquid, the pressure in the header is computed using 
the perfect gas law as follows: 


1 T = (M nc + M v + Mi). + (M nc + M v + M ] ) st 


R Mt T 


r ' 144 M W V 0 
Where : 

V 0 = volume in header, ft J 

If the pressure P calculated is greater than the input down- 
stream pressure, then iteration ceases and the mass of gas 
leavinq the header is computed to make the mass of gas remain- 
1ng V give the downstream pressure. To do this the follow, ng 
relations are used: 
p dx 

M e = Mj (1 - p-) 

m S t = m T * ir- 
on the other hand, if the pressure P calculated is less than 
the downstream pressure, the iteration is repeated with pres- 
sure P used in HBALNC instead of the downstream pressure. The 
iteration repeats until the pressure tried equals the pressure 
calculated using the perfect gas law. When this occurs, a 
the gas stays in the header. Therefore: 

M e = 0.0 

M S T = My 
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3.5.1 


The Headers (Continued) 

The case for the gas being a mixture of non-condensables and 
vapor was just described. The following describes the analyses 
for the case where the gas is all vapor. First, analyses are 
performed assuming all the steam will stay in the header be- 
cause the downstream pressure is too great. In this case, 
temperature is iterated upon until the energy equation is 
bal anced . 


Energy In + (Energy Stored)i = (Energy Stored )2 


Hi - Mi hgj 
h ST,1 = M ST,1 


h 9,ST 


R (Tst + 460) 


(1) 

( 2 ) 


m t = M STj1 
HsT, 2 = m T 


+ Mi 

[hg,TRY - 


R 1 

; t tf 

}Y + 460) 

i 

(m w : 

>v J 


h ST, 2 = H ST,1 + H i 


(3) 

(4) 


Temperature Typy is iterated upon until equation (4) is satis- 
fied. At this point, the pressure in the header is computed 
assuming all the vapor stays in the header and using the new 
temperature of stored vapor in the header . 


t ST, 2 = t TRY 

M t R (Tst,2 + 46 °) 

P = H4 (M w ) y V 0 

If this pressure is indeed less than the downstream pressure, 
then all the steam does stay in the header. If this pressure 
is greater than the downstream pressure, then the energy bal- 
ance must be repeated with an exiting enthalpy term added. 

Energy In + (Energy Stored)! = Energy Out + (Energy Stored )2 


_ 144 P d (M w ) u V 0 
R (Typy + 460) 

M e = My - M$T,2 

HST,2 = MST.2 


R (Tjry + 460) 
h 9,TRY - (Mw)”j 


H e = M e hgjRY 

H e + HsT, 2 = H i + H S T,1 


(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 
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Entering masses 


Figure 8 

SAWD Bed Energy Balance 
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Figure 9 

control Volume Around Bed Segments 
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3 . 5.1 The Headers (Continued) 

Tomnprature Ttrv is iterated upon in equations (5) through (9) 
2ST* r 5Stl2fW i« satisfied. At that po nt the^exlting^ 
gas temperature and temperature of ga . ■] 

set equal to T T RY- The Pressure in the header is set equ 

the downstream pressure. 

3 # 5 # 2 The Resin Plus Void Plus Ca nister^ 

e .»ir.nr. 

Figure 8. 

A nf oas M- enters the void volume at temperature Tj during a 

Ji^e step t! \h mass of gas is assumed to mix with the ex, ting 
aas M in the void volume which is at a temperature Bi* , 

matching 1 numbered segments is a mass of water and carbo ■ 

"Bj numbered S&X'l £st ttW&Y 

mixture and the bed Tb 2 . The new mass of gas in the void volume 
Mu and its temperature^? are used in the perfect gas law to ca 
rui^tP the Dressure in the segment void volume if all the gas were 

dtoxid^and 0 water Z .TiM ^ plus 

foam M r are all at the final temperature 

The energy balance equation begins with the uniform state uniform 
flow form of the first law of thermodynamics as presented in Re 
erence 6 by Van Wylen and Sonntag: 


Q 


( u ; 


V 2 



t .i.) 

h 9c / 

= W + 

^M e (h e 

V 2 
+ T- 

2gc 

2 

- + Z 2 
c 


Z M i 

( n 2 4 
\ U1 + V 

• Zl -3- 
1 9c 


+ lc 


r) 


) 


control volume around the 


I old volljme^ segment'and' the' corresponding numbered resin segment as 

shown in Figure 9. The assumptions in applying this first law are. 
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The Resin Plus Void Plus Canister (Continued) 

1) The control volume remains constant relative to the 
coordinate system. 

21 The state of the mass crossing the control surface is 
’ constant with time and uniform over the various areas of 
the control surface where flow occurs. 

31 The state of the mass within the control volume may change 
’ w??h time: but, at any instant of time, the state is uni- 
form over the entire control volume (or over several iden- 
Jlflable regions that make up the entire control volume). 

For the SAWD control volume, as presently modelled. 


W = 0.0 
0 = 0.0 


Vo = Vi = 0.0 

Z 2 = Zi = Z e = Zi = 0.0 


A, so, it is assumed that changes in kinetic energy are ne 

to^i temizef each^of°the r consti tuents^and 1 then^ rearranged te a fL 
suUable for the computer. The first law can be expanded to: 


0 = (M Vj e h v ,e - M v,i h v,i) + ( M nc,e h nc,e - M nc,i h nc,i) + 

(H c g ,e heg.e - Mcg.i heg.i) * <Mv,2 "v.2 - N,,1 u v ,i) * 

(M nc ,2 Unc ,2 - M nc,l “nc,l) * ( M cg,2 “eg, 2 - N cg,l M cg,l u cg,l) * 
{M w ,2 u„, 2 - Mw,i u„l) * (M C 1,2 del ,2 - M cl ,1 “cl.l) + 

(M r ,2 u r , 2 - H r ,i u r ,l) + <"t,2 u t ,2 - «t,l “t,l + 


(Mf ,2 u f ,2 “ M f,l u f,l^ 

The following mass balance relations are invoked for the non 
condensables : 


M nc,i + M nc,l " M nc,e + M nc,2 
AM nC = M nc>2 - M nc,l 
M n c,e = M nc,i - ^ M nc 

The following mass balance relations are invoked for water 


M v ,i + M v ,l + Mw.l = M v ,e + M v ,2 + M w ,2 
A M w = M w > 2 - M w ,l 

AM V = M v> 2 “ M v ,l 

M v , e = M v ,i - AM W - AM V 
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Th» Resin Pl us Vojd Plus Can ister (Continued) 

~ : r^T^ions are invoked for carbon dioxtde 

The following mass balanc 

AM C 1 = M C 1 ,2 - Mcl ,1 
AM C g = M cg , 2 - M cg ,1 

M o = Mca i ■ ^M c l ’ AMcg 

C9, • ss balance relations are invoked for the resin, 

Th e following mass balance 

tank, and foam: 

M r , 2 = M r,l = M r 
Mt,2 * M t,l = M t 
Mf,2 = M f,l = Mf 

Invoking these «,ass relations yields: 

0 = M v ,i (bv.e-^v.il^nc.i One,. - h nc ,i * 

<W (W - *cg,i> * M ».l (u »- 2 • UV ’ 1) + 

M nc ,l (u„c,2 - Rnc.l) + M cg,l ( u cg,2 - “«>1 

M U ,1 (u u , 2 -u„,l)tHcl,l (901,2-^1,1)* 

M r (u r ,2 - u r ,l) + M t (dt,2 - ut.l) (9f,2- Uf . 

AM„ (hv,e - 9W,2) - A H cl (M,e-“cl,2)- 
AM V (b v ,e - u v ,2) - ^"C9 (M.e ' u cg,2> ‘ 

AM nC (h nc ,e " u nc , 2 ) 

Now specific heats and temperatures are introduced where: 

T v , e = T nc ,e = T cg,e = T e 
T v ,i = T nc ,i = T cg»i = Ti 
T v ,2 = T nc ,2 = T cg,2 = T2 
T v ,1 = T nc,l = T cg>l = Tl 
T w 2 = T C 1 ,2 - T r ,2 ‘ T ». 2 * Tf - 2 = Tb * 2 
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ThP Resin Plus Void Plus Ca nister (Continued) 
Substituting temperatures and specific heats yields. 

0 = M Vji c P)V (T e - T i) + M nc,i c P,nc ( T e - T i ) + 
(T e - Ti) + M Vj l c v , v ( t 2 - T l) + 


M 


'eg »i C P > c 9 

M nC( l c v ,nc (T2 - T l) + M C9>1 Cv ’ c 9 ( T2 ' Tl) + 

M w ,i c w (T b , 2 -T b ,l) +M cl,l c cl (Tb ,2 
M r c r (Tb,2 * T b,l) + M t c t (T b,2 * T b. l} + 

M f c f (Tb,2 - Tb.l) - ( h v,e - u w,2) - 
AM cl (h cg>e - u cl , 2 ) * aM v ( h v,e " u v,2) " 

AM cg (h cg , e - u cg , 2 ) - AM nc ( h nc,e - u nc,2) 

j f thp exitinq conditions from the 
^ 1n the vo,d se9ment - 

Therefore: 


T 2 = T e u v ,e = u v,2 

h v ,e = hv,2 u nc,e = u nc,2 

h nc ,e = hnc,2 u cg,e' u cg,2 

heg ,e = h cg ,2 

Substituting : 

0 = ( M v,i c p,v + M nc,i c p,nc + M cg,i c p,cg) ( T e ^ 

(M v> i c v , v + M nc> i c v , nC + M cg,l c v,cg) ( T e “ T l } + 

(M w ,l c w + M c1 ,1 c cl ) (T b ,2 - T b,l) + 

(M r c r + M t c t + Mf Cf) (Tb,2 * T b,l) 

AM W (h v , e - u W ,2) * (h ^9» e ' Ucl ’ 2) " 

AM V (h v , e - u Vl e) - AM cg ( h cg,e - u cg,e) - 

AM nc (hnc ,e “ u nc,e) 


But : 

144 P v RT 
h " u = J = M w 
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The Resin Plus Void Pl us Canister (Continued) 

0 = (M v ,i c p ,v + M nc,i c p,nc + M cg,i c P,cg) ( T e ' T i } + 

(M v ,1 c v,v + M nc,l c v ,nc + M cg,l c v,cg) ( T e 
(Mw,i c w + Mcl.l cel) ( T b,2 - T b,l) + 

(M r c r + M t ct + MfCf) (Tb,2 - T b,l) - 

aM w (h V) e - u w,2) " AM cl ( h cg,e ' u cl ,2) ' 

(AM v /Mw v + A M cg /Mw c g + AM nc /Mw nc ) R T e 

But: 

b v ,e - u w ,2 = h fg + Ah «a 

Where : 

A h wa = heat released daring water absorption, Btu/lbm 

It is assamed that Ah w » is insignificant compared to h fg and 
is therefore ignored. Therefore. 

b v ,e - U W,2 = bfg 
A1 so : 

s? sk '■ ” 

val ue of 1270 Btu/lbm. See Reference 8. 

Substituting these two relations yields: 

0 = (M v j c pjV + M nc> i cp, nc + M cg,i c P,cg) ( T e Tl ^ 

(M v> i c v ,y + M nc ,l c v ,nc + M cg,l c v,cg) ( T e * T i ) + 

(M w ,l c w + M c i,i c c i) (T b ,2 “ T b,l) + 

(M r Cp + M t c t + M f Cf) Ob, 2 " T b,l) 

AM W hfg - AM C 1 a h cr - 

(AM v /Mw v + A M C g/Mw C g + AM nc /Mw nc ) R T e 
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Thp Resin Plus Void Plus Caniste r (Continued) 

Lastly, it is assumed that: 

T e = ^b,2 
T 1 = T b,l 

0 = (M vJ c P jV + M nc>i c p>nc + M cg ,i c p ,cg ) ( T e ‘ T i 
(M v>1 c v>v + M nc> i c v ,nc + M cg,l c v,cg) ( T e T i ^ 

(M w ,i c w + M c i f i c c i) (Te - T i) + 

(M r c r + M t c t + Mf c f ) (T e “ T i ) * 

AM W hfg - ^M c ] Ah cr - 

(AM V /Mw v + AMcg/MWcg + AM nc /Mw nc ) R T e 

In this equation, there are six unknowns: 

T = gas exit and final bed temperature, °F 

aMw - change in the mass of absorbed water on the bed, Ibm 

aMr , . change in the mass of absorbed carbon dioxide on the 
bed , 1 bm 

A „ v - change in the mass of vapor in the void volume segment, 
lbm 

A „ nc - change in the mass of non-condensables in the void 
volume segment, lbm 

A M, a = change in the mass of carbon dioxide gas in the void 
y volume segment, lbm 

The procedure - 

5T5.VjrsrJC xrrsr - 

ss siV.- 5 -rs aa sss zxsszax 

is determined from the following equations. 

M c i = r M r j^( p C 02^ i 1 ' ^ P C02^eq] At/3600 
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The Resin Plus Void Plus Canister (Continued) 

The absorption/desorption rate r is calculated from the 
following equation: 

r = r 0 (1 - 0.9 L H 2 o) 

Where: 

r 0 = baseline absorption/desorption rate = 1.5 pph/ 
lb IR45/psid 

l_ H 20 = bed water loading, lb of M per lb of dry amine. 

The equilibrium partial pressure of carbon dioxide (P c0 2) e q 
is calculated from the following equations: 

< P C02>eq ' V C02 P T 

Where Ypn? is a function of equivalent carbon dioxide loading. 
This function is derived from data in Reference 7 and Hamilton 
Standard's own experimental data. 

The equivalent carbon dioxide loading is determined from the 
actual carbon dioxide loading and two adjustment factors as 
fol 1 ows : 

(L ) = ^02 

1 C02'eq WF ACT * TFACT 

WFACT represents the effect of water loading on the carbon 
dioxide loading. The effect is determined from Hamilton 
Standard test results and shows that actual CO 2 loading Lqq 2 
increases with water loading and reaches a maximum around 20 
to 30% water loading. This agrees with conclusions reached 
by P. Mignon in Reference 7 where he states: 

"Below 10 percent, water deficiency reduces the possibility 
of reaction between amine and dissolved CO 2 . Above 30 
percent, the water layer is too important and limits the 
CO 2 diffusion toward aminated molecules." 

TFACT represents the effect of temperature on the carbon 
dioxide loading. The effect is determined from Hamilton 
Standard test results and shows that actual CO 2 loading J-C02 
decreases significantly with increases in temperature. This 
agrees with data by P. Mignon in Reference 7 and data presented 
by A. Boehm in Reference 9. 

Therefore, to determine the mass of carbon dioxide absorbed or 
desorbed, the parameters that need to be known are: 
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3.5.2 


The Resin Plus Void Plus C anister (Continued) 

T e = gas exit and final bed temperature, °F 

Mp = mass of dry resin, lbm 

(P ) = partial pressure of carbon dioxide in the total 

^C02'il !j! he inlet an d state 1 gases. 


At 

LH20 

PT 

TO 


= time step, seconds 
= water loading 
= total pressure, psia 

= baseline absorption desorption rate, pph/lb 
dry amine/psid 

Therefore, during a given 


Of these, the only unknown is T e - 
time step: 

AM c i = f (Te) 

<Ve2 * leg’ll ‘ A " cl 

The amoont of carbon ^ iWnXiiT 

Of carbon dioxide aval lab le desorbed , s limited to 

gas; while the amount of carbon dioxl °| s ,ding in the resin. 

the amount of carbon dioxide ^ g(j Qn the car bon dioxide 

Lastly, the following minimum P a before any flow 

left in the void segment after absorption 
He ralrnlated to leave the void: 


144 (^cojjea ^002 v ° se g 

^cg^e2,min R (T e + 460) 

After solving the a ^? ve . n^iVthe bed AM^the following 

equation^ar^sol ved 1 to 1 detennine cVnli in the mass of 

absorbed water on the bed aMw. 

AMw = (M v ) n - (M v ) e2 

The total mass of vapor (M y ) e2 to be ^ft in void segment 
and to exit the void segment is determined fr om^the foil o wl 

(M v ) e2 = Mw v [n nc ,i + n n c,l + ( M c gl e 2 /fMwc9 ^ P T - p H20,eq 
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th. Resin Plus Void Plus Canister (Continued) 

. . / m \ -j c known from the previous calcu- 

In th,S ** ,2’ n diSSldl The only unknown is the partial 

lations on carbon dioxide. y . } preS sure is 

pressure of the water, P H20 ,eq- lh P 

determined from: 

(Weq = Psat (Te) 1f (P H20 ) il^ P5at (Te) 

(Weq - RKeo k Tsat (Te) " < Wil i ^ ^ 

( P H20>eq ' (P H20>il 1f RHeq X Psat ^ < Ph2 ° 11 S *\ 

Where the equilibrium relative humidity RH e , Is determined from 

SSiTty 

ibrium testing. The model *??!£!? ‘ „Yh the bed water and 
of the exiting gas is in e ?“’ ’ '5"™ "uming that the 

carbon dioxide loadings. In addition by ass ^ pressure 
loadings are constant over a time step, P vo lume 

can MSiSS TS^Tof'Z, ?he final bed temper. 

ature. 

The total mass of vapor of 

and to exit the void s ®9 me ceament plus that originally in 

vapor (Mv)il entering the 9 nriainallv in the dry resin 

the void segment plus the mass M w ,l originally in y 

itself. In other words: 

(M v ) e2 < Hv,1 + M v,l + M « al 

Like AM cl , the only unknown to determine AM W is the temper- 
ature T e * 

A My, = f (T e ) 

M nc,2 = ( M nc) e 2 

M v,2 = ( M v^ e 2 

M cg,2 = ( M cg) e 2 

AMpc = ^nc ,2 " ^nc ,1 
AM V = M v>2 - M V> 1 
AM c g = M c g,2 - M cg ,1 
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The Re sin Plus Void P lu^janister (Continued) 


3 I I I ' *"* *** _ , „ _ 

ad < p then a pressure collapse 
However , if P T < p d a " d P } h - S ^ not a possible situataon " , 
situation has occurred. Th logic is used to adjust 

reality. In this event, t e 0 '« of ca rbon dioxide 

the changes in the m ?sses ( A cl i y too much carbon dioxide 
and water on the than was physically possible 

or water was absorbed onto tn void volume. 

k°™ IS.’lSfltrVS rXSV downstream pressures. 

144 P min Vo 
^min “ R (T e + 4°^ 

(M ) = M v ,2 Main - ( M nc>e 2 /M "nc ' M “v 

W V this point in the 5rtemined.° 

water and carbon dioxnde on ti^bedjuve , p ^ vo1d segme „t 

Therefore, the total determine d and are: 

plus those which exit are 

(H nc ) e2 - M nC ,2 + M nc,e = M nc,l + M nc,i 

(M v ) e 2 = M v ,2 + M v,e = Mv>1 + Mv,i W 

,» i s m 9 + M cg , e = M cg,l + M 09.i " AMcl 
M °9 e2 c9, ■ h is the split between the mass of 

What remains to be determined^ and the mass which 

exi^^Then'wllenThis split 1 ^Xted^o^lli^he energy 

^quati^^ nC the n spl^t C ietween e the^mass^remain^ng^^n^the ped 

to^maintain tSeY-i^reater than the downstream pressure. 
The Pressure that -ould ex ^ segment is give n by: 

(M v ) e2 , and the (M cg ) e2 remai 

Mr R (T e + 46Q ) 

P T = 144 Mw T Vo 

Where: 

W T = (M nc ) e2 t («cg>e2 + lb ” 

Mwj = molecular weight of My 
Vo = segment void volume, ft 

P T hypothetical pressure in void volume, psia 
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The Resin Plus Void Ph ^anister (Continued) 

If the pressure P T > Pd, then the splits are given by: 

Mnc,e - ("nc) e2 x 

M v ,e - ( M v) e 2 x (1 ‘ Pd/Pl) 

M C g,e = ( M cg) e 2 x ^ ' Pd/Pj ^ 

M n c ,2 = ( M nc ) e 2 x Pd/Pj 

M v ,2 » (M e 2 X Pd/PT 

M c g ,2 = ( M cg) e 2 X P(J//PT 
AM nC = M nc ,2 " M nc,l 
AM V = M Vj 2 " M v,l 
AM c g = ^cg ,2 " M cg ,1 

p < but P T > Pm, then all the mass stays in 
If the pressure Pj < p d Dul - n _ u 
the void segment and: 


M nc ,e = M v,e = M cg,e 


am 


0.0 


- M v ,2 - M V,1 
AM W = (M v ) n - ( M v ) e 2 

If only carbon dioxide were absorbed, then: 

(Mcg)e 2 - "eg, 2 - M. ' («"c) e2 /"“nc - («») e2 /"“vl "“9 
AM C g = M cg> 2 - M cg ,1 
AM c i = (Mcg ) i l " ^ c 9^e2 

If both carton dioxide and water were absorbed, then: 

N m in - ( M nc) e 2 /Mw " c 


F = 


TmTT^Tm^" (McgTTi/^cg 


(M t ) e2 

- F * (H v ) n 

( M cg>e2 

= F * (M cg ) n 

AM V 

= M Vj2 - «v,l 

AM c g 

= M c g, 2 " M cg ,1 

am c i 

■ ("cghl - (M cg> 

A Mw 

= ( M v ) i i - ( M v)e2 
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3.5.2 The Resin Plus Void Plus Canister (Continued) 

Now the values for aM v , aM w , AM C g, AM C ], AMnc» an< ^ ^e are 
substituted into the energy equation. If the sum is less than 
0.1, the energy equation is assumed to be balanced. If not, 
temperature is iterated upon using the standard Newton-Raphson 
technique. However, if the next trial temperature is greater 
than or equal to the saturation temperature for the downstream 
pressure, the temperature to try will be permitted to approach 
the saturation temperature in 25 increments from the initial 
temperature to the saturation temperature. 


3.5.3 


T e 


Tj (25 - I) + T S at * \ 
25 


Thus when I = 25, T e = T sa t- Once the temperature T e reaches 
T sa t» then temperature is no longer iterated upon but instead 
the mass of vapor (M ) 2 remaining in the void segment plus 
exiting is iterated up8n. 


The Temperature Sensor 


The temperature sensor is at the exit of the SAWD bed and has a 
characteristic time constant. The standard differential equation 
for the response of a temperature sensor to a fluid temperature is: 


dT 

dt 



Where : 

= temperature sensor time constant, sec 

T a = ambient temperature of fluid surrounding the 
temperature, F 

T = temperature sensor temperature, F 

In difference form the differential equation becomes: 
T a - T 1 

T 2 - Ti = 3 - — x At 


Or : 

T a " T 1 

T 2 = Ti + ~ ■ At 

This equation is used to determine the temperature as indicated 
by the temperature sensor. Currently, the time constant has a 
value of 240 seconds. 
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3.6 


Tanks 

The TANKG subroutine is used for ^^^ z l^ k the A W compl et^descri ption 
well as the carbon dioxide accu R ^^ al^he analytical 

of the subroutine is given in Reference 3 for an 

options. 

For the water supply tank, the following options are selected: 


1 ) 


The outlet flow during changes 

S°« during steady state 

to : 


m £ 


= 2.55 lbtn/hr 


2 ) 

3 ) 

4 ) 


for steady state operation. 

No transient heat balance is performed; instead, the steady 
state heat balance equations are used. 

Initial mass = 168.3 lb 
Initial pressure = 30.0 psia 
Initial temperature = 70 F 

i) no provision for heat rejection/addition through use of secon- 
dary flow as heat transport medium exists. 

luring transient operation the flow out : of tank Is determined^ 

SS.’®! Mil be described in 
Section 3.12.3 for SAWD Control. 

For the carbon dioxide accumulator tank, the following options are 
sel ected : 

11 ["he'infet n^desj.?^^"^^ transient, ‘the 

GP0LY1 coding will be executed. 

" SS?^Swm=3s 

for steady state operation. 
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Tanks (Continued) 

3) No transient heat balance is performed; instead, the steady 
state heat balance equations are used. 

4) Initial volume, pressure, and temperature inputs are used to 
calculate the initial mass. 

Initial volume = 2.0 ft^ 

Initial pressure = 30 psia 
Initial temperature = 70°F 

5) No provision for heat rejection/addition through use of secon- 
dary flow as heat transport medium exists. 

During transient operation, the flow out the tank is regulated by 
orifices, valves, and pressure regulators as shown in Figure 2. If 
the tank pressure ever exceeds 30 psia, the back pressure regulator 
opens to vent the excess pressure to the overboard vent. Below 30 
psia, the flow from the carbon dioxide tank can be directed to the 
overboard vent or to a carbon dioxide reduction process. If dir- 
ected overboard, the solenoid valve opens at 30 psia and closes at 
24 psia. The pressure in the tank will therefore oscillate up and 
down between 24 and 30 psia. If directed to a carbon dioxide re- 
duction process, the corresponding solenoid valve is opened at 24 
psia and the flow is regulated to exit at 0.275 lbm/hr by the down- 
stream pressure regulator and orifice combination. Below 24 psia, 
no flow exits the tank and tank pressure is allowed to build. This 
exit flow logic is represented mathematically by the following 
equations : 


For P > 30 psia: 


m 


e 


. 144 Mw Vo (P - 30) 3600 

m- + = i L 

1 R (T + 460) At 


Where: 


Mw = 44.01 for carbon dioxide 
m e = exiting mass flow, lbm/hr 
m, = inlet mass flow, lbm/hr 

The above equation is designed to permit sufficient exit 
flow during the time step At to return the tank pressure 
to 30 psia. 

For : 


P < 24 psia : 
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3.6 Tanks (Continued) 

144 Mw Vo (30 - P) 3600 
m e = mi * R (T + 4bif) At 


If: 

m e < 0.0, m e = 0.0 

The above equations are desi qnedtoperm' t0 
there is no inlet flow and also “ t T4 ensure from ever 

MTdWpliW "■ ,n one t1me ste ?; t o T r h1s restrict1on 
simulates the 30 psia back pressure regulator. 

For pressures between 24 and 30 psia and flow directed to a 
rarhon dioxide reduction process. 


m e = 0.275 lbm/hr 

For pressures between 
overboard vent: 


24 and 30 psia and flow directed to an 


m e = 0.2316 x P 

This relation is 
of 0.065 and the 


derived from a solenoid valve 
following flow equation: 


flow coefficient 


SCFM = 13.61 P 


C v 


\J SG x (T + 460) 


Where: 

p = pressure 

C v = valve flow coefficient = 0.065 

SG = specific gravity = 1.529 

T = temperature = 70°F 

SCFM = flow at standard temperature (70°F) and 
pressure (14.7 psia), cfm 


_ _ SCFM „ M 

m p — — “ X l 

e 386.73 
Substituting all 


m e = 0.2316 P 


x 60 x P 

the appropriate values yields. 
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Tanks (Continued) 

Now both expressions for the exit flow for the tank Pressure 

between 24 and 30 psia are restricted so * s one 

pressure fron exceeding 30 psia or falling below 24 psia in one 

time step. 

The minimum exit flow for pressures between 24 and 30 psia is 
given by: 

144 Mw Vo (30 - P) 3600 
1545 (T + 460) At 


m e,min m i 


The maximum exit flow for pressures between 24 and 30 psia is 
given by: 

144 Mw Vo (P - 24) 3600 
m e,max = m i + 1545 (T + 460) At 


Pumps 

The PUMP subroutine is used to analyze component 28 the water supply 
pump A complete description of the subroutine is given in Refer- 
ence 3 for all the analytical options. 

For the water supply pump, the only calculation performed is to add 
the input pump heat to the fluid stream. No head rise or °t"er Pump 
performance calculations are performed. The exiting temperatur 

given by: 


T e = 


Ti + 


Q 


m y c p 


The 


Where 6 i s the heat input to the fluid stream by the pump, 
pump heat addition used in ESCM is 10.0 watts. 

Steam Generator 

The SMGEN subroutine is used to analyze the component 29 steam ” 
generator. A complete description of the subroutine is given in 
Reference 3 for all the analytical options. 

For the steam generator in ESCM, heat is supplied by electrical 
heating, and the subroutine computes the amount of e ]® c i r J? a ; e e 9y 
reou red to vaporize the entering fluid and to superheat it to a 
user specified amount. The amount of electrical power required is 

calculated as fol lows: 
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3>8 stem Generator (Continued) 

Where: 

m.: = incoming mass f1ow ’ 

h 1 = exiting enthalpy, Btu/lbm 

h® = inlet enthalpy, Btu/lbm 

h-j = h-) (Ti ) 

Whe re : 


V» M C I ^ • 

(Ti ) - enthalpy of saturated liquid at temperature 
1 1 i; T i# Btu/lbm 

h e = bg (T e ) 

Where: 

h (T ) , enthalpy of saturated vapor at temperature 
9 6 T e » Btu/lbm 

T e ■ T sat ( p ) + AT spbt 
Where : 

T sa t ( p ) = saturation temperature at pressure P, °F 

at nht = desired degrees of superheat, F 

• h hpat to electrical power, the following 
To convert the required heat to electr 

relation is used: 


n / O A 1 7 


3.9 


fitters 

— subroutine H used to 

SPLIT subroutine is used as follows. 


Component 

No. 


3 Cabin secondary 
heat exchanger. 

9 Humidity control 
water separator 


Component Descr iption 
flow splits to SAWn and humidity control 


heat exchanger primary flow splits to 
and water separator bypass. 
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3.9 


Splitters (Continued) 


Component 

No. 


rnm ponent Desc ription 


14 

17 

22 

24 

26 

30 


, .. sAWD bed #1 and to 
Cabin secondary flow spins to SAW 

I J » o 

, ... c A uo fan and to C0 2 accura- 
SAWO bed *1 exit flow spins to SAW 

ulator or bed #2. 

i-+c to SAWD fan and to CO 2 accum- 
SAW bed n exit flow spins to SAW 

ulator or bed #1. 

f1 to CO? accumulator or 

cab i p* and* to° t he^ op po s 1 te° n umbe red 

Stea. 9 enerator or SAWbed « or .exit flow splits to 
SAWD bed #1 and #2 ml 

SAWbed*! or « exit flow splits to cabin and C0 2 acco»- 
u 1 a to r • 


u 1 a to r • 2 

M e c^» 

... rafin f or component number 3 is aexen 
sultinVsecondary°and tStal flows: 


m se c = ( m e) 19 + ^ m e^24p 

m T = (n e )i2 + ( m e)lO s + 


(m e )24 


= m sec /m-|- 


SR 

Where : 

( % ) 19 = SAWD fan exit flow, lbm/hr 


(■I,” . splitter * 2 4 exit priory flow, lb./br 

-I ^ r* fan PX 


\ a = spl ~ - 

;, 2 p . nonidit, control beat exchanger fan exit flow,lb./br 

(m e )io s * " ater se P arator " ater fl0W ’ 

= split ratio 
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3.9 Splitters (Continued) 

This adjustment of the split -tie prior to solving compo -^number 

th^cabi ^humidit^contro^heat ex^anger and to the SAWD carbon 
dioxide removal subsystem. 

input data set. The values used in ESCM are. 


SR V = 0.03 
SR] “ 1 .00 

SRo2 = °- 03 

SRn2 = 0.03 


SRc02 " 0.03 
SRjC = 0.03 

SRsF#l = °- 03 


3.10 Mixer 

The GASMIX subroutine is used to mix Tse^ 

i n°Fi gS?2 e 3 h tnd 3 ?tb1 e lTthe GASMIX subroutine is used as follows: 


Component 

No. 


Component Description 


11 Flows from the water separator and the water separator 
bypass line mix. 

13 Flows from the humidity control heat exchanger fan and the 
SAWD desorption flow sensor mix. 

15 Flows from the cabin and the steam generator or bed #2 mix 
before entering SAWD bed #1. 

18 Flows from SAWD beds #1 and #2 mix before going to the 
fan . 

20 Flows from the cabin and the steam generator or bed #1 mix 
before entering SAWD bed #2. 

23 Oesorption flows from SAWD beds #1 and K mix before going 
to the opposite bed, back to the cabin, or to the C0 2 

accumul ator . 

25 Steam generator flow and flow from SAWD beds #1 or #2 mix 
before going to the inlet of one of the SAWD beds. 

31 Flow from the cabin mixes with returning flow from the SAWD 
before entering the humidity control heat exchange . 
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3.10 


3.11 


Mi xer (Continued) D , 

A complete description of . the D ^assumed ^^"temperature . 
ence 3. An adiabatic mixing process a mixed stream are 

dew no nt, pressure, and constituent f1ows J ts 15 and 20, the 

7 related . In all instances except 0 r secondary 

pressure is taken as the ^^^from the source which has a non- 
source pressure or as the P -, wayS take the pressure of th ® 

zer0 to maintain pressure in the SAWD bed after 

desorption S flow ends and before absorption begins. 


All Purpos e Components 
The ALTCOM subroutine is used to handle the ^^exchangers and to 
EH^^^r:r U t;^r^ption of At™ is 
given in Reference 3. 

When used as secondary or 

humidity control heat | J* * a repository for the cooling 

executed. The componen * *JJ to the heat exchanger. The humidity 

flow temperature , and pre nressure are constant at 

cSoling fl0 ": d te J e t te F„ P r ESCh( the values are: 

input values specified by the user. 

Flow = 950 lbm/hr 
Temperature = 45°F 
Pressure = 50 psia 

The sensible heat exchanger cooling f l ~s^tiSSi?a 2 %^HSever^ 
STeSSSt^rlSK and pressure are constant for ESCM at 

and 50 psia respectively. no 2 hereby the heat 

The water separator ' s ^ a '^^ n “ S {emperatu?e is calculated. Prior 

ALTCOM , the following as per- 

formed in GP0LY1. . 

1} Any entrained liquid from the condensing heat exc an 
removed from the primary flow. 

2) The separator heat load is added to the flow stream. 

q = 3.84 x 3.413 - 13.1 Btu/hr 

. macc flow of vapor, and mass 
In ALTCOM, the assuming an adiabatic mixing. 

See W a°description ^subroutine HBALNC in Reference 3. 

i f rt is that all the entrained liquid entering the 
water S separator 3 is by the water separator. 
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3.12 Control 

Three items are controlled in the ARS subgroup analyzed by ESCM. 

They are the cabin temperature, the cabin oxygen and total pres- 
sures, and the SAWD carbon dioxide removal subsystem. The control 
logic is written in subroutine GP0LY1. Here will be presented the 
control logic for each item. 

3.12.1 Two Gas Controller 

The two qas controller for ESCM controls the entering flow of oxygen 
and nitrogen into the cabin to maintain the desired oxygen pressure 
from 3 09 to 3.23 psia. The controller admits nitrogen to bring the 
pressure up to 14.813 psia only after the oxygen partial Pressure 
has been brought above 3.09 psia. The logic is shown in Table Z 
while the opening and closing flow schedules for oxygen and nitrogen 
are shown in Figures 10 and 11, respectively. These curves are 
based on July 1980 Pressure Control Subsystem testing at Johnson 
Space Center. 

3.12.2 Cabin Temperature Control 

The cabin temperature is controlled by the cabin sensible heat 
exchanger. To remove more heat, the cooling flow to the heat 
exchanger is increased; and conversely to remove less heat, the 
cooling flow is decreased. The calculation of the proper cooling 
flow is done in GP0LY1 using a Newton-Raphson iteration technique 
described in Reference 3 in the subroutine ESTIM. The basic equa- 
tion is: 

m 2 “ m l IT T \ 

m = mi + A — (T \ - T se tl 

T 1 - >2 


Where : 

m = new cooling flow, lbm/hr 
mi = flow at state 1, lbm/hr 
Ti = temperature at state 1, °F 
m 2 = flow at state 2, lbm/hr 
T 2 = temperature at state 2, °F 
A = iteration constant = 0.025 

The minimum flow is limited to 950 lbm/hr while the maximum 
flow is limited to 1900 lbm/hr. 
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When flow calculated by closing curve = 0.0, valve closes. 
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CABIN TOTAL PRESSURE (PSIA) 


Figure 10 

Oxygen Opening And Closing Flow Schedules 
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3.12.3 SAWD Control 


jr vn u ~ 

The sequence of operation is '^^^^q^orbs'whi 1 e'bed 0 *! is on 
first full cycle. «***"?{**£ Tr £ea» flows through the 
standby. During standby, ... bed ^ 2 is directed first to 

bed. Initial flow from J he . d ^°^ J the CO? accumulator with 
the cabin to prevent con ^ 1 "f J SlVnfas detected by the flow 
oxygen. When the flow r Vq at this point, is directed to 

sensor, the flow, | lm ost total y 002 t acc ^ ulat0 r until 24 psia 
the accumulator. Pressure bul1d * it ted to leave at a 

is reached. At this P^ essU j e ’ ti nrocess If pressure builds 

sssw. s;<M bri° 

pressures K-rW. SS P-ents any bach 

flow from the accumulator to the beds. 

Oesorbing in bed ,2 continues until Us «« ^and 

180°F , which indicates most C0 2 b J s Doint steam flow ceases 
flow consists mostly of steam. allow the pressure in bed #2 

and appropriate valves are sw' f 0 e e conds. Then Abe fan starts 
to bleed slowly to bed #1 for id d then to the humidity 

pulling cabin air through bed iZtobe* ^ through bed #2 it 

control heat exchanger ’nlet. A * s V ie^ and thereby heat s bed 

pushes the steam in bed #2 f five m i nu tes or until the 

#1. This energy transfer ve first, 

bed #1 exit temperature reaches 130 r , wmcne 

After energy transfer 1 S h C ?T P lh1 Te * bed G #2 begins absorption. Valve 
ected to bed #1 to desorb J be d #: 2 and air exiting 

positions are changed to . i heat exchanger inlet. Desorb- 

the bed flows to the humidity contro^ back t0 t he cabin via 

ing continues in « " th ^ eachcs 0 .01 cfm. Then, the desorb 
the flow sensor until t , , The Drocess continues as 

flow is directed to the accumulator. and bed #2 absorb- 

described previously, but with be # .. ^ bed #2 continues for 

ing. For only this startup the average cabin relative 

a time which is determined directly from^ne ^ ^g R< The t1me 

humidity measured during the p dete f m1ned from the average cabin 
for all other absorption y rhe previous absorption cycle. 

MI then^al ternat^ through" thi ^sequence as directed by the 

SAWD controller. 

One aspect that needs ^^cycirto^l^dependlSg on 

The desorption steam . ? tbe cab i n during the past absorp- 

S, th the Ca followi ng rel ation : 

m . 0.0396 (49 + 0.32867 RHj - 0.008671 RH { 2 ) 
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4.0 


SAWD Control (Continued) 

Where: 

RH % = average cabin relative humidity in percent 
nomenclature 


Document 

Symbol 

C v 

H 

J 

L 

M 

MR 

Mw 

N 

P 

Q 

0 

QL 

QLMIN 

QS 

QT 

R 

RH 

SCFM 


Description 

Valve flow coefficient 

Total enthalpy, Btu 

Joule's constant = 778.26 ft-lbf/Btu 

Bed loading, lb/lb-dry amine 

Mass of substance, lbm 

Per man metabolic heat load, Btu/man-hr 

Molecular weight, lbm/ lb-mole 

Number of people in crew 

Pressure, psia 

Heat transferred to system, Btu 
Rate of heat transfer, Btu/hr 

Latent heat added to gas stream per man, Btu/man-hr 

Minimum latent heat added to the gas stream per man 
Btu/man-hr 

Sensible heat added to the gas stream per man, 
Btu/man-hr 

Total gas stream heat load from the crew, Btu/hr 
Universal gas constant = 1545 ft-lbf/mole-R 
Relative humidity 

Volumetric flow at standard temperature (70°F) and 
pressure (14.7 psia), cfm 
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4.0 


NOMENCLATURE (Continued) 


Document 

Symbol 

SG 

SR 

T 

TCAB 

U 

V 

Vo 

W 

WC02 

WH20 

W02 

Z 

cfm 

C P 

c v 

h 

m 


Description 

Specific gravity with respect to air at 70°F and 14.7 
psi a 

Split ratio in SPLIT subroutine 

Temperature, °F 

Cabin air temperature, °F 

Total internal energy, Btu 

Velocity, ft/sec 

Volume, ft ^ 

Work done by system, Btu 

Total crew carbon dioxide generation rate, lbm/hr 

Crew water vapor generation rate, lbm/hr 

Crew oxygen gas usage, lbm/hr 

Elevation of system, ft 

Volumetric flow, cubic feet per minute 

Constant pressure specific heat, Btu/lbm 

Constant volume specific heat, Btu/lbm 

Specific enthalpy, Btu/lbm 

Mass flow, lbm/hr 


n 

r 


Lb-mol es 

Reaction rate constant for absorption of carbon diox- 
ide pph of CO? per lb of dry resin per psi difference 
between pressureof CO? in void segment and equilibrium 
pressure in void segment. 

Base reaction rate constant at zero water loading, pph 
of CO, per lb of dry resin per psi difference between 
pressure of CO? in void segment and equilibrium pres- 
sure in void segment. 
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NOMENCLATURE (Continued) 


Document 

Symbol 


Descri ption 


Specific internal energy, Btu/lbm 
Power, watts 


Mole fraction 


Greek Letters: 


Subscripts: 


Change in mass of carbon dioxide gas in segment void 
volume, lbm 

Change in mass of carbon dioxide in resin, lbm 

Change in mass of non-condensables in segment void 
volume, lbm 

Change in mass of vapor in segment void volume, lbm 
Change in mass of water stored in resin, lbm 
Pressure difference, psid 
Temperature difference, °F 

Heat released during carbon dioxide absorption, Btu/lbm 

Heat released during water absorption, Btu/lbm 

Time step, seconds 

Heat exchanger effectiveness 

Fan efficiency 

Density, Ibm/ft^ 

Time constant, seconds 
Absolute humidity 


Carbon dioxide 


Water 
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NOMENCLATURE (Continued) 


Document 

Symbol 

Subscripts : 

02 

N2 

SF#1 

ST 

T 

TC 

TRY 

a 

b 

eg 

cl 


Description 


Oxygen 

Nitrogen 

Special flow #1 

stored 

Total 

Trace contaminants 

Trial value during iteration 

Ambient or air 

Bed 

Carbon dioxide gas 

Carbon dioxide'stored in resin 


d 

db 

dp 

e 

eq 

e2 

g 

i 

il 

1 

min 


Downstream 

Dry bulb 
Dew point 
Exit 

Equilibrium 

Sum of exit plus state 2 conditions 1 

Gas or saturated vapor 

Inlet 

Sum of inlet plus state 1 conditions 
Entrained liquid or saturated liquid 
Mi nimum 


n void segment 


in void segment 


57 




UNITED 

TECHNOLOGIES 

[K)MfflOO=ir®K] 


SVHSER 9504 


4.0 NOMENCLATURE (Continued) 


5.0 


Document 

Symbol Description 

nc Non-condensables 

r Resin 

sat Saturation condition 

spht Superheat 

t Tank or metal canister around SAWD resin 

u Upstream 

v Water vapor 

w Water stored in resin 

1 State in control volume at start of time step 

2 State in control volume at end of time step 
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